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We discuss the future prospects of heavy neutrino searches at next generation lepton colliders. 
In particular, we focus on the planned electron-positron colliders, operating in two different beam 
modes, namely, e^e“ and e“e“. In the e^e“ beam mode, we consider various production and 
decay modes of the heavy neutrino {N), and find that the hnal state with e + 2j + $, arising from 
the e^e“ —>■ Nv production mode, is the most promising channel. However, since this mode is 
insensitive to the Majorana nature of the heavy neutrinos, we also study a new production channel 
e^e“ —>■ Ne'^W^, which leads to a same-sign dilepton plus four jet hnal state, thus directly probing 
the lepton number violation in e"'"e“ colliders. In the e“e“ beam mode, we study the prospects of 
the lepton number violating process of e“e“ —>■ W~W~, mediated by a heavy Majorana neutrino. 
We use both cut-based and multivariate analysis techniques to make a realistic calculation of the 
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We hnd that with the cut-based analysis, the light-heavy neutrino mixing parameter |14Afp can be 
probed down to ~ at 95% C.L. for the heavy neutrino mass up to 400 GeV or so at y/s = 500 
GeV with 100 fb~^ of integrated luminosity. For smaller mixing values, we show that a multivariate 
analysis can improve the signal signihcance by up to an order of magnitude. These limits will be at 
least an order of magnitude better than the current best limits from electroweak precision data, as 
well as the projected limits from y/s = 14 TeV LHG. 
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I. INTRODUCTION 


The observation of neutrino oscillations in solar, atmospheric, reactor and accelerator neutrino data 
(for a review, see e.g. [1]) has unequivocally established that at least two of the three active neutrinos have 
a non-zero mass and that individual lepton flavor is violated. This provides a conclusive experimental 
evidence for the existence of some new physics beyond the Standard Model (SM), and therefore, a precise 
understanding of the neutrino mass mechanism is an important step in unraveling the nature of new 
physics. 

The minimal renormalizable extension of the SM to explain the observed smallness of neutrino masses 
is the neutrino portal, defined by the Lagrangian 

-C = h,^L,^N„ + ^MN^^N^Np+Il.c., (1) 

where = {vi, £)~[ (with £ = e,p,T) and $ are the SU{2)l lepton and Higgs doublets respectively, 
<I) = 1(72$* ((72 being the second Pauli matrix) is the hypercharge conjugate of <f>, = N^C~^ [C 

being the charge conjugation matrix) denotes the charge conjugate field and Na are SM gauge-singlet 
neutral fermions, usually known as the sterile neutrinos, since they can talk to the SM sector only via 
their mixing with the active neutrinos. In Eq. (1), hia are the dimensionless complex Yukawa couplings 
which, after the electroweak symmetry breaking, yield a Dirac mass matrix Md = hv, v being the Higgs 
vacuum expectation value. In addition, Eq. (1) allows a Majorana mass term Mjv, since the sterile 
neutrinos carry no SM gauge charge. The existence of these new interactions in Eq. (1) leads to the 
following neutrino mass matrix in the flavor basis {vf ,Na\: 


M, 


( 0 Md\ 
\Ml Mn) 


( 2 ) 


The active-sterile neutrino mixing is parametrized by ^ ~ In our subsequent discussion, we 

often denote the elements of this matrix as Vm (see Section H for details). For ||^|j = \/Tr(^t^) <C 1, 
the light neutrino masses and mixing are given by the diagonalization of the effective mass matrix 

M, ~ -MdM^^mI , (3) 

whereas the heavy neutrino masses are of order M^. This is known as the type-I seesaw mechanism [2-6]. 

From the above discussion, it is evident that there are two key aspects of the seesaw mechanism that 
can be probed experimentally: the Majorana mass Mjv of the sterile neutrinos and their mixing Vin 
with the active ones. The Majorana nature of both active and sterile neutrinos can in principle be 
probed via the lepton number violating (LNV) low-energy process of neutrinoless double beta decay 
(Oiy/S/S). However, an observation of Oi/PP does not necessarily probe the active-sterile mixing since 
these mixing effects may not always give the dominant contribution over purely left-or right-handed 
contributions to the process; for a review, see e.g. [7]. Alternatively, a non-negligible value of 

VfN can be indirectly inferred from non-unitarity of the light neutrino mixing matrix [8, 9] as well as in 
observables for lepton flavor violation (LFV) [10, 11], lepton non-universality and electroweak precision 
tests [12-16]. However, these low-energy observables by themselves do not prove the Majorana nature 
of neutrinos since models with pseudo-Dirac neutrinos can also give rise to large non-unitarity and LFV 
effects [17-22]. 

On the other hand, the direct collider searches for heavy neutrinos at the energy frontier can provide 
a simultaneous probe of both the aspects of the seesaw mechanism, i.e. the Majorana nature of the 
heavy neutrinos N and their mixing with the active neutrinos, provided the seesaw scale is accessible to 
the available center of mass energy; for a review, see e.g. [23].^ In particular, electroweak-scale heavy 
neutrinos can be produced on-shell at colliders with a cross section depending on the strength of their 
mixing with the active neutrinos in the minimal setup, while their subsequent decay pattern is governed 
by whether they are Majorana or Dirac particles. At a hadron collider, the ‘smoking gun’ collider signal 
for heavy Majorana neutrinos is the same-sign dilepton Hnal state with two jets and no missing transverse 
energy: pp —)■ W* —)■ —>■ £^£^jj [41-46]. Using this channel, both CMS and ATLAS experiments 

at the LHC have set direct limits on the light-heavy neutrino mixing |V^Ar|^ (with £ = e, p,) for heavy 
neutrino masses between 100 and 500 GeV [47-50]. With the run-H phase of the LHC with more energy 


^ For some examples of low-scale seesaw models with relatively large active-sterile neutrino mixing, which lead to observable 
effects at colliders as discussed here, see [19, 24-40]. 
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and higher luminosity, and including the infrared enhancement effects due to t-channel photon-mediated 
processes [51], these limits could in principle be extended for up to 1 TeV or so. 

For electroweak-scale Dirac neutrinos, the same-sign dilepton signal is suppressed, and the ‘smoking 
gun’ signal at the LHC is the trilepton final state: pp —> W* —>■ N£^ —>• [52-57]. Analyzing 

the trilepton data from -v/s = 8 TeV LHC [58], similar limits on |V>ArP (with £ = e, u) have been derived 
for Mn up to 500 GeV [57]. 

In light of vigorous discussions (see e.g. [59]) on the possibility of a future lepton collider, such as 
ILC [60], FCC-ee [61], CLIC [62], CEPC [63] and muon collider [64], it is worthwhile examining the 
prospects of heavy neutrino searches in the clean, almost background-free environment of a lepton collider. 
We recall that for heavy neutrino masses below the Z-boson threshold, using their possible production 
in the Z-boson decay Z —)• viN or Z —)• Nv^ [65], and its subsequent neutral current (NC) or charged 
current (CC) decays, 95% C.L. limits on [V^Ar]^ < 10“^ were obtained by L3 [66] and DELPHI [67] 
collaborations from a reanalysis of the LEP data. A future high-luminosity Z-factory, such as the FCC- 
ee, will significantly enhance the sensitivity by several orders of magnitude [68]. Proposed fixed-target 
experiments such as SHiP [69] can also dramatically improve the limits on IV^Ar]^, though only in the 
low-mass range between 0.5 and 10 GeV or so. 

For heavy neutrino masses above the Z-boson threshold, a direct search was performed in e“'"e“ 
annihilation at LEP [70, 71], where a single heavy neutrino can be produced via its mixing with active 
neutrinos: e+e“ —)■ Nv [72-77], and can subsequently decay via GG, NC or Yukawa interactions with 
the SM W, Z or Higgs (ff) boson respectively: N —)• £~W^ (£+lV“), v{v)Z, v{v)H. Concentrating on 
the decay channel N —)■ elF with W —)■ jets, which would lead to a single isolated electron plus hadronic 
jets and missing energy, the L3 collaboration put a 95% C.L. upper limit on the mixing parameter [VeAfP 
in a heavy neutrino mass range between 80 and 205 GeV [71]. Similar limits were derived recently [16] 
using the ALEPH data on —)■ W~W~^ —>■ v£~£'^v [78]. However, these searches were kinematically 

limited by the maximum center-of-mass energy ^/s = 208 GeV at LEP2. It is expected that future lepton 
colliders can significantly improve the sensitivity and extend the limits to higher heavy neutrino masses 
for both Dirac and Majorana cases. This is explicitly demonstrated here for the most promising channel 
e+e” —)■ Nv —>■ £jj$'^ by performing a Monte Garlo simulation at ^/s = 350 and 500 GeV ILG, including 
realistic detector effects and background estimates. We have used both cut-based and multivariate 
analysis (MVA) techniques to derive the projected sensitivity limits for a generic ILC detector [79]. We 
find that a light-heavy mixing of [VeAr]^ down to ~ 10“^ can be probed at 95% C.L. for the heavy 
neutrino mass up to 400 GeV or so at ^/s = 500 GeV ILG with 100 fb“^ integrated luminosity. Although 
for concreteness, we have performed our simulations for the ILG benchmark parameters, this study is 
also relevant in the context of other future lepton and photon colliders [80]. 

Note that neither the existing LEP search channels e+e” —> Nv —)■ £jj$ nor —)■ W~W~^ have 

LNV final states, and hence, it is not evident from these channels whether the heavy neutrinos are 
Majorana particles, as predicted by the usual seesaw mechanism. In principle, the Majorana or Dirac 
nature of the heavy neutrinos could be determined to some extent from the e“'"e“ —>■ Nv —>■ £jji$ mode 
using the distribution of the polar angle between the produced neutrino and the incoming electron [81]. 
Here, we propose an alternative clean way to determine the Majorana nature using a new LNV signal at 
e“''e“ colliders within the minimal seesaw: 

e+e" ^ Ne^W^ £^e^ + Aj , (4) 

with same-sign dilepton plus four hadronic jets. Noting that the SM background for this process is 
extremely small, we perform a detailed detector-level simulation for the signal sensitivity and find that 
the detection of 10 signal events for a mixing parameter [14^1 = 0.04 would require an integrated 
luminosity of at least 700 fb~^, which is within reach of next generation lepton colliders. 

For completeness, we also present the heavy neutrino production cross sections for various other 
subdominant processes, such as the Higgsstrahlung process e+e“ —> ZH —)• ZNvi{vi) and the dominantly 
photon-mediated processes e+e“ —>■ N£^Vf.{ve). We could also have the processes —> NW^ and 

e ^7 —>■ N£^Ve{ve), with a real photon coming from the electron or positron beam. These processes are 
of direct relevance for an electron-photon collider [82]. 

Finally, taking advantage of the fact that it is possible to switch the beam configuration of an e“'"e“ 
machine from e+e” to an e“e“ mode [83], one can also look for an LNV signal in the process e“e“ —)■ 
W~W~ —)■ 4j mediated by a t-channel Majorana neutrino [84-88]. We make a detailed detect or-level 
analysis of this LNV signal and find that due to negligible SM background for this process, it offers 


^ For a lepton collider, since the initial four-momenta of the particles are known, we can measure the full missing energy 
f and not just the transverse component, as in the case of a hadron collider. 
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an interesting alternative way to probe LNV in the electron sector, even though the cross section is 
suppressed by |T4Ar|^. 

The plan of the paper is as follows: in Section II, we discuss various possible production modes of the 
heavy neutrino at a lepton collider. In Section III, we do a detailed event analysis using the cut-based 
techniques for the signal and background involving the final states with e + 2j -\- $ and -I- 4j at 

the e+e“ collider and the LNV 4j final state at the e~e~ collider. In Section IV, we present an MVA to 
further enhance the signal sensitivity for the most promising channel with e + 2j + final state. Our 
conclusions are given in Section V. 


II. HEAVY NEUTRINO PRODUCTION AT LEPTON COLLIDER 


In the minimal seesaw scenario, where the heavy neutrinos do not carry any gauge charge, the only 
way for them to communicate with the SM sector is through their mixing with the active neutrinos, 
as governed by the seesaw Lagrangian (I). To parametrize the light-heavy neutrino mixing, we first 
diagonalize the full neutrino mass matrix given in Eq. (2) by a unitary mixing matrix: 

V~^MuV = diag(TOi, Mj) , (5) 


where mi (with i = 1,2,3) and Mj (with j = 4, 5,...) are respectively the light and heavy neutrino mass 
eigenvalues. The unitary matrix V has an exact representation in terms of a dimensionless matrix ^ 
(which depends on Md and Mn): [37, 89] 




( 6 ) 


where 7/, 7/' are the unitary matrices diagonalizing the individual light and heavy neutrino mass matrices 
My and M^v respectively. 

Now using Eqs. (5) and (6), the light neutrino flavor eigenstates can be related to the mass 
eigenstates 9i and Nj as follows: 


ve = 




U' 






lAfiijNn = UiiVi Vi-iNi 


( 7 ) 


where the first term on the right-hand side (RHS) measures the non-unitarity of the Pontecorvo-Maki- 
Nakagawa-Sakata mixing matrix, and the second term determines the size of the light-heavy neutrino 
mixing in CC and NC interactions involving neutrinos. In the charged-lepton mass diagonal basis, the 
CC interaction relevant for the production and decay of heavy neutrinos is given by 

-Ccc = -^W-IrPL^e + ^-c. = + +H.C. , (8) 


where Pf, = (I — 7^)/2 is the left-chirality projection operator. Similarly, the NC interaction is given by 


—Cnc = 


^ ^ ^ . Z„ 


2 cos 9. 


2 cos 


+ iv^v),jNa^PLN, 


+ +H.C.} 


( 9 ) 


where 9^, is the weak mixing angle. For our subsequent discussion, we use a model-independent phe¬ 
nomenological approach, parametrized by a single heavy neutrino mass scale Mjv, assuming that any 
other heavy neutrinos present in the system are sufficiently heavy and do not affect our analysis. Hence¬ 
forth, we denote the light-heavy mixing parameter simply as without explicitly writing the heavy 
neutrino flavor index. This enables us to derive generic collider exclusion/sensitivity limits in the 
(Mjv, Vm) parameter space, without referring to any particular low-scale seesaw model-building aspects. 

From (8) and (9), we see that the cross section of heavy neutrino production in e“'"e“ collision will be 
proportional to the mixing parameter jV^Arj”, where n = 2 or 4, depending on the process considered (see 
below). Similarly, for the two-body decays of the heavy neutrino to the SM final states ^ Zv{v) 

and Hv{v) (if kinematically allowed), the corresponding partial decay widths are also proportional to 
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FIG. 1. The branching ratios (BR) of the heavy Majorana neutrino decay modes N —>■ Zu{u), Hv{u) 

above the Z-boson threshold in the minimal seesaw model. 


\Vm?- 
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4 
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^ Ml 


. ^Ml\ 

l + 2T7f 

Ml 


( 10 ) 


( 11 ) 


( 12 ) 


Note that the total decay width for a heavy Majorana neutrino will be twice the sum of the partial 
decay widths shown in Eqs. (10)-(12), after taking into account the charge conjugate processes. The 
corresponding branching ratios (BR) are shown in Figure 1. For larger values of Mff, the branching 
ratios follow BR(.^1F) : : BK{Hiy) ~ 2 : 1 : 1. Note that BR(7V —)■ £~W~^) in the heavy 

Dirac neutrino case is twice as large compared to the heavy Majorana neutrino case, which has equal 
probability to decay into either £~W'^ or t^W~. 

In what follows, we consider various heavy neutrino production channels at a lepton collider and 
numerically compute their production cross sections for given heavy neutrino mass and center of mass 
energy.^ For our subsequent collider analysis, we have implemented the relevant Lagrangian terms 
as shown in Fqs. (8) and (9) for heavy Majorana neutrino interactions in FeynRules2.0 [90], and 
generated the Universal FeynRules Output (UFO) [91] model files. These UFO model files have been 
used in the the Monte-Carlo (MC) event generator MadGraphS [92] to generate the parton-level 
cross sections shown in this section, as well as the signal and background events for the analysis in 
Section III. The following basic trigger cuts on lepton transverse momentum (pp), pseudo-rapidity (1) 
and radial distance {AR^ were used to identify the leptons in the final state: 


Pt 


> 20 GeV, \l\ < 2.5, Ai?" > 0.4 


(13) 


For our numerical simulation, we consider two benchmark values for the center of mass energy ^/s = 350 
and 500 GeV. Also, we focus on heavy neutrino masses only above the Z-threshold. For Mm < Mz, 
the existing LFP limits from Z-decay are quite stringent [66, 67[ and one needs a high-luminosity Z- 
factory [68, 93[ to significantly increase the sensitivity. 


A. 






A single heavy neutrino of any flavor can be produced through an s-channel Z-exchange in e' e 
collisions, as shown in Figure 2 (a).'^ The corresponding cross section will be proportional to [V^ArP- For 


® For general analytic expressions for the production cross section, see [74]. 

^ If A is Majorana, both and are allowed. However, for brevity, we denote both as simply ui in the following. 
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FIG. 2. (a) The s-channel and (b) i-channel Feynman diagrams for the process e''"e —>■ Nui. 




(a) (b) 


FIG. 3. Normalized heavy neutrino production cross section for the channel e''"e“ —>■ uiN for = 350 and 500 
GeV. The left panel corresponds to £ = e (both s and t channels), whereas the right panel corresponds to i = fj,,T 
(only the s channel). 


i = e, there is an additional t-channel FF-exchange process, as shown in Figure 2 (b). The differential 
cross section for this process can be found in [73, 75]. The numerical values of the total production 
cross sections are shown in Figure 3 for two different values of y/s = 350 and 500 GeV. The left panel of 
Figure 3 corresponds to the case £ = e and the right panel is for the case i = fi,T. 

The enhanced cross section in the electron channel is due to the additional t-channel contribution. 
Similarly, for the heavy neutrino decay, N —>■ eW has the largest branching ratio, as evident from 
Figure 1. Thus, the channel e“''e“ —)■ Nve, followed by iV —)■ eW and W —)■ jets, has the largest signal 
cross section, and hence, can act as the most promising channel for the discovery of heavy neutrino 
mixing with electron neutrinos at a lepton collider. As a matter of fact, the existing direct search limit 
from LEP [71] comes from this channel. 

Due to the relative smallness of the s-channel contribution, this production mode does not provide a 
very promising signal for probing IVfArP with £ = Also note that the s-channel contribution to the 
cross section decreases with increasing y/s, except when the heavy neutrino mass is close to the kinematic 
threshold, where the phase space suppression becomes more dominant. This is evident from Figure 3 
(right panel). Using a suitable beam polarization might increase the signal sensitivity for £ = fi,T [94], 
but this is of limited practical interest for a small VeN- 

We note here that since the hadronic activity at an e“'"e“ collider is very limited, the background can 
be easily controlled with the help of simple kinematic cuts. Therefore, the other decay channels of the 
heavy neutrino, i.e. N —>■ Zvi, Hvi, can also be used as complementary search channels, even though 
the corresponding BR are about a factor of two smaller than the N ^ W£ mode. This will be illustrated 
in Section IV. 


B. e+e- ^ 

The signal discussed in Section IIA is insensitive to the Majorana nature of the heavy neutrino, 
and hence, does not probe LNV at a lepton collider. Here we present a new production mechanism 
e+e” —> N£'^W^ that can probe the Majorana nature of the heavy neutrinos at an e+e” collider. This 
is due to the fact that for a Majorana neutrino, both £'^W^ decay modes are allowed; hence, we get 
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FIG. 4. Some sample Feynman diagrams for the process e^e“ —>■ N£~W^, as discussed in Section IIB. The 
diagrams of type (a) can probe light-heavy neutrino mixing in all flavors, whereas the diagrams of type (b) and 
(c) can probe only the electron-sector mixing. 




(a) (b) 

FIG. 5. Heavy neutrino production cross section for the channel e''"e“ —>■ for y/s = 350 and 500 GeV, 

and with \VfN\ = 0.04. The left panel is for mixing with electrons, whereas the right panel is for mixing with 
muon and tau sectors. 


a LNV same-sign dilepton signal, along with two Ws. Thus, the dominant signal in this case will be 
+ Aj which probes LNV. In addition, for £ e, one can in principle also probe the LFV at future 
colliders. 

This process gets contribution from several diagrams, and three typical diagrams are shown in Figure 4. 
For the electron sector, there is a dominant contribution from the t-channel photon diagrams [cf. Figure 4 
(b)], thus leading to an infrared enhancement effect.® This enhancement effect is absent in muon and 
tau sectors for an e“''e“ collider. The total production cross section for this process is given in Figure 5. 
The left panel shows the cross section for electron-sector mixing and the right panel is for € = /i,r. Note 
that this process involves diagrams with a heavy neutrino as as intermediate state [cf. Figure 4 (c)], and 
therefore, the total cross section for this process does not simply scale as jV^ArPi since the total heavy 
neutrino decay width in the propagator explicitly depends on [FfArp. The cross section values shown 
in Figure 5 were obtained assuming a typical mixing parameter jV^ivI = 0.04. A detailed detector-level 
simulation is done in Section III to derive the sensitivity for this process at future lepton colliders. 

A related process to consider is —>■ NW^, where the real photon comes from one of the electron 
or positron beams. This process can be calculated using the Weizsacker-Williams equivalent photon 
approximation (EPA) [96, 97]. The results are shown in Figure 6. Again in this case, one can have a 
same-sign dilepton signal, when N —?> and the other W with the same sign as the lepton in the 

final state decays leptonically. 


C. e+e- ^ 

Similar to the case discussed in Section IIB, the process e+e” —)■ also contains infrared 

enhanced t-channel contributions mediated by virtual photons, as illustrated in Figure 7(a). The cross 
section for this process is given in Figure 8. Similar to the situation in Section IIB, this process involves 


® This is similar to the infrared enhancement effects in hadron collisions: pp —>■ Ni^jj at the LHC [51, 57, 95]. 








Mn (GeV) 


FIG. 6. Heavy neutrino production cross section for the channel e^7 —>■ NW^, with the real photon coming 
from one of the electron or positron beams, for y/s = 350 and 500 GeV, and with IVejvl = 0.04. 




FIG. 7. Sample Feynman diagrams for the process e''"e —>■ as discussed in Section IIC. 




(a) (b) 

FIG. 8. Heavy neutrino production cross section for the channel e^e“ —>■ for y/s = 350 and 500 GeV, 

and with mixing \Vin\ = 0.04. The left panel corresponds to i = e and the right panel is for i = jj,, whereas (! 
could be any of the three lepton flavors. 


Feynman diagrams with heavy neutrino propagators [cf. Figure 7(b)], and hence, the cross section does 
not simply scale as 

A related process to consider here is —)■ where the real photon comes from one of the 

electron or positron beams. We calculate the cross section for this process using EPA and the results are 
shown in Figure 9. 


D. e+e- ^ ZH ZNvt, 

This is the Higgsstrahlung process, followed by the decay of Higgs to as shown in Figure 10. 
However, this process is only effective for values below the Higgs mass, as evident from the cross 
section values shown in Figure 11. This process can directly probe the heavy neutrino mixing with all 
neutrino flavors. Note that for the heavy neutrino mass in the vicinity of the Higgs mass, the Higgs 
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FIG. 9. Heavy neutrino production cross section for the channel e *7 —>■ Ni^ve, with the real photon coming 
from one of the electron or positron beams, for y/s = 350 and 500 GeV, and with II^jvl = 0.04. 



FIG. 10. Feynman diagram for the process e^e —>■ ZH —>■ ZNvi. 



FIG. 11. Normalized heavy neutrino production cross section for the channel e''"e —^ ZNve for = 350 and 

500 GeV. 


decay width will get significantly modified due to the new decay channels available, and using the Higgs 
signal strength data, one can also derive indirect constraints on the mixing parameter [98, 99]. 


E. e e —>■ W W 



FIG. 12. Feynman diagram for the t-channel N exchange. 
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FIG. 13. Normalized cross section of the process e e —>■ W W for y/s = 350 and 500 GeV. 


Experimentally, it is possible to switch from e“'“e“ mode to the e~e~ mode within the same experi¬ 
mental setup [83]. In this case, the Majorana nature of the heavy neutrino can be tested by observing 
the LNV process e~e~ —?■ W~W~, as shown in Figure 12. This process is mediated by a t-channel JV 
exchange,® and is free from SM background, except in the case when two additional leptons are produced 
and lost in the beam pipe. This can be considered as the lepton-collider analogue of the Qvfif} process. 
Moreover, for a linear collider like the ILC, it is also possible to use polarized electron beams. By using 
two left-handed electron beams, the cross-section for the e~e~ —)■ W~W~ channel will increase by a 
factor of four as compared to two unpolarized beams. On the other hand, for a circular collider like the 
FCC it is not possible to retain the polarization of the colliding beams due to the strong magnetic field 
used to bend the beam tracks. 

The cross section for the channel shown in Figure 12 is proportional to and the normalized 

cross section is shown in Figure 13 for ^/s = 350 and 500 GeV. The produced Ws can decay either 
hadronically or leptonically, thus generating the final states with Aj, 2£~ + $ or 2j + l~ + ^. The cross 
section values after the selection cuts and the signal sensitivity will be discussed in Section III. 


III. CUT-BASED ANALYSIS 

We have used the MC event generator MadGraph5 [92] to generate all signal and background events. 
The subsequent parton showering and hadronization have been done using PythiaO [100]. We have used 
Delphes3 [101] for simplihed detector simulation relevant for a generic e“''e“ collider. In this section, 
we present a cut-based analysis (CBA) to discriminate signal from the background and to derive the 
sensitivity limits for heavy neutrinos. In Section IV, we have employed a more sophisticated MVA for 
better signal-to-background discrimination, which leads to a better signihcance. 

For our subsequent collider analysis, we apply the following basic trigger cuts, in addition to those in 
Eq. (13), to identify leptons and jets in the final state: 

> 20 GeV, |77^ | < 2.5, > 0.4, > 0.4. (14) 

Jets are clustered using the anti-fer algorithm [102] implemented in FastJet [103] with clustering pa¬ 
rameter R = 0.4. Then we apply further selection cuts for each signal discussed below, depending on its 
kinematics, in order to isolate it from the SM background, if present. 

In the following, we consider various hnal states corresponding to the different heavy neutrino produc¬ 
tion modes, as discussed in Section II. Since the cross sections involving the mixing with the muon and 
tau sector are much smaller compared to the electron sector, we only consider the latter to derive the 
discovery prospects of the heavy neutrino at a future lepton collider. 

A. e“*“e“ —>■ iVzZg 

The dominant SM background for this final state comes from the e+e” —>■ W~^W~ channel, followed 
by the leptonic decay of one W and the hadronic decay of the other. Therefore, it is expected that 


® A similar diagram exists with a t-channel light Majorana neutrino exchange. However, its contribution to the cross 
section is extremely small, suppressed by the light neutrino mass. 
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Mn (GeV) 

(a) 

FIG. 14. Significance defined in Eq. (16) as functions 
and (b) -s/s = 500 GeV, with C = 100 fb“^. 
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Mff for different mixing values at (a) y/s = 350 GeV 


the invariant mass of the electron and the neutrino, viz. M{e,^) peaks around the W-mass for the 
background, while for the signal, one does not expect such a behavior. On the other hand, the invariant 
mass distribution of the isolated electron and the two hardest tagged jets, M(e, ji, j 2 ), shows a peak near 
Mjv for the signal, but not for the background. These two distributions for the signal and background 
are shown later in Figure 18 in the context of a MVA. Therefore, M(e, and M(e, ji, J 2 ) distributions 
can be used to separate the signal from the background in order to get a better sensitivity. Based on this 
observation, we demand the following selection cuts specific to this signal: exactly one isolated electron 
and at least two tagged jets satisfying the basic trigger cuts, no photon with pT > 20 GeV and 

\M{e,ji,j 2 )- Mn\ < 40 GeV, \M{e,$)-Mw\ > 20 GeV. (15) 


Mn 

(GeV) 

Signal (fb) 

Background (fb) 

Afs 

TG 

SC 

TC 

SC 

V-^s +Ab 

250 

21.94 

18.96 

638.25 

68.99 

20.22 

350 

16.33 

13.48 

638.25 

46.11 

17.46 


TABLE I. Signal ( 0 - 5 ) and background (ug) cross sections after the trigger cuts (TC) and the selection cuts (SG) 
as defined in the text for different Mn with IVejvl = 0.04 at yT = 500 GeV. The significance is computed using 
100 fb“^ integrated luminosity. 

After imposing these selection criteria and including the detector effects, we compute the signal and 
background cross sections for two illustrative values of the heavy neutrino mass Mjy = 250, 350 GeV at 
y/s = 500 GeV and for the mixing parameter |14Ar| = 0.04. The results are shown in Table I. The signal 
significance 


" = vmrm ■ 

where JVs and JVb are the number of signal and background events respectively, has been computed with 
an integrated luminosity of 100 fb“^. 

In Figure 14, we show the significance n defined by Eq. (16) as a function of the heavy neutrino mass 
Mat, for different mixing values |14Af| = 0.04,0.03 and 0.02 at ^/s = 350 and 500 GeV for C = 100 
fb“^. We see that mixing values up to \VeN \ = 0.02 can be probed with a significance of > 5a for heavy 
neutrino mass % 350 GeV at ^/s = 500 GeV. 

In order to make a comparative study of the projected sensitivity obtained here with the existing 
and future limits from other contemporary experiments, we show the 95% C.L. projected ILG limits 
on iVeATp as a function of Mn in Figure 15. We have considered two configurations with i/s = 350 
(green curves) and 500 GeV (magenta curves) and with integrated luminosities C = 100 fb“^ (solid) and 
500 fb“^ (dashed). For larger Mat, the cross section decreases, as evident from Figure 3, thus leading 
to a weaker limit on |t4Afp- For comparison, we have also shown the current 95% C.L. direct exclusion 
limits from LEP 2 (red, shaded region) [71] and ^/s = 8 TeV LHC data (blue, shaded region) [50], as 
well as a projected 95% C.L. limit from ^/s = 14 TeV LHC with 300 fb“^ integrated luminosity (blue, 
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(GeV) 

FIG. 15. The projected 95% C.L. ILC sensitivity of the light-heavy neutrino mixing parameter |I4jvP as a 
function of the heavy neutrino mass Mn for = 350 (green) and 500 GeV (magenta), and with 100 fb“^ 
(solid) and 500 fb“^ (dashed) integrated luminosity. Also shown are the current 90% G.L. indirect exclusion 
limit from electroweak precision data (EWPD; brown, dotted) [16], 95% G.L. direct exclusion limits from LEP 
2 (red, shaded) [71] and ^/s = 8 TeV LHG data (blue, shaded) [50], as well as a projected 95% C.L. limit from 
= 14 TeV LHC with 300 fb~^ integrated luminosity (blue, dashed) [23]. 


dashed curve) [23]. The best current limit is the indirect limit at 90% C.L. (brown, dotted line) derived 
from an analysis of the EWPD [16].^ As can be seen from Figure 15, a linear collider can significantly 
improve the heavy neutrino sensitivity up to mass values close to its kinematic threshold. Note that here 
we have not shown the constraints from 0;//3/3 (see e.g. [33, 36, 104-107]), since these limits are highly 
model dependent and can be significantly weakened in certain cases due to cancellation between different 
contributions [108]. 

As far as the other decay modes of the heavy neutrino are concerned, it can go to Zv and Hv final 
states, followed hy Z ^ H ^ bb and W ^ iv decay modes, thus giving rise to jj + 

bb + $ final states. For illustration, we show in Table II the signal cross sections of these 
final states for two benchmark values of Mm at ^/s = 500 GeV. Note that the + $ final state has 
two contributions coming from the leptonic decay of either W or Z, as shown by the last two columns 
of Table II. Due to the relatively clean environment at a lepton collider, these processes can still lead to 
a non-negligible number of events with 500 fb~^ luminosity. Note that for the channel involving Higgs 
decay, we have taken into account the correct Higgs decay width at next-to-next-to-leading order (NNLO): 
V{H —>■ 66)nnlo = 0.0023 GeV [109], whereas the corresponding width obtained from MadGraph5 at 
leading order is 0.0053 GeV. 


^ For a detailed discussion of the present bounds on light-heavy neutrino mixing, see [23]. 
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|y„jv| = 0.04; 500 GeV 


FIG. 16. (a) cross sections for the process e''"e“ —>■ Ne^W^ —>■ e^e*+4j after trigger cuts and including detector 
efficiency at ^/s = 350, 500 GeV and (b) the required luminosity to observe 10 signal events with ^/s = 500 GeV. 
Here we have taken the mixing \VeN\ = 0.04. 


Mn (GeV) 

2j + $ (fb) 

2h + $ (fb) 

2£ + $ {Z) (fb) 

21+ $ (IT) (fb) 

250 

16.41 

9.03 

1.641 

11.05 

400 

7.415 

5.20 

0.7417 

4.926 


TABLE II. Parton level signal cross sections for the subdominant decay modes of the heavy neutrino produced in 
the e+e“ —>■ Nvi process after trigger and selection cuts for illustrative cases of heavy neutrino masses Mm = 250 
and 400 GeV and for mixing parameter \ VtN\ = 0.04. 


B. e“*“e —>■ Ne^W^ + 4j 

In this subsection, we consider the LNV process of e+e” —>■ Ne^W^ —>■ + 4j, as discussed in 

Section IIB. We demand in the final state, two electrons (or two positrons) and at least four tagged 
jets satisfying the basic trigger cuts given in Eqs. (13) and (14). The signal cross section after imposing 
the cuts and including the detector efficiency is given in Figure 16 (a). The SM background for this 
process, mostly coming from 41 + 4j final states with two leptons lost in the beam pipe, is negligible 
< 0(10“® fb). Hence, we do not perform a detailed background simulation for this process. It is evident 
from Figure 16(a) that the signal cross section is also small (< 10“^ fb) and drops sharply with increasing 
Mat. This makes it a challenging task to observe a reasonable number of signal events with 100 fb“^ 
integrated luminosity. In Figure 16 (b), we show the required integrated luminosity for observing 10 
signal events. We find that C > 700 fb“^ is required for this purpose at a ^/s = 500 GeV ILC. This 
is nevertheless within reach of future lepton colliders [80]. Although this channel has a smaller cross 
section, as compared to the e + 2j + channel discussed in Section III A, an observation of even a 
few LNV events in this channel will concretely establish the Majorana nature of heavy neutrinos, with 
profound implications on our understanding of the origin of neutrino mass. 


C. e-e-^W-W-^4j 

In e“'"e“ colliders, the positron beam can be readily switched to an electron beam, thus converting 
to an e“e“ machine. This freedom allows us to study LNV processes such as e“e“ —)■ W~W~ —)■ 4j, 
as discussed in Section IIF. At the partonic level, there are no SM backgrounds for this LNV process. 
However, in realistic situations, different other processes can mimic the same final state as the signal. For 
instance, we can have 4j + 4j + $+ leptons or 4j+ leptons, where the additional leptons are lost in 
the beam pipe. The SM processes which contribute to the background are listed in Ref. [87]. We revisit 
these processes here and implicitly assume that we cannot reconstruct the charge of the gauge bosons 
from where the jets originated. The dominant background channels under consideration are : (PI) 
e“e“ -4 W~W~ vivi, (P2) e“e“ -4 W~ Zvie~ , (P3) e“e“ —>■ ZZe~e~ and (P4) e“e“ -4 W~W'^ e~e~ . 
All the gauge bosons decay hadronically so that we have four jets at the parton level. We shower the 
sample events and perform the detector simulation by demanding a minimum of 20 GeV, no leptons 
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FIG. 17. The significance as a function of Mn for the e e —>■ 4j channel for an integrated luminosity of 500 
fb"^ at ys = 350,500 GeV. 


and photons which pass the trigger cuts and at least 4 tagged jets satisfying the trigger cuts given in 
Eq. (14). The background cross sections before and after the selection cuts are shown in Table III. For 
•y/s = 500 GeV and |V),iv| = 0.04, we can achieve > 3cr significance for the full range of the heavy neutrino 
mass considered here, as shown in Figure 17. However, given the smallness of the cross sections, the 
significance is somewhat lower at y/s = 350 GeV. We comment here that a more involved search strategy 
using optimized cuts or a multivariate analysis might improve this significance. Also note that, in our 
analysis, we have used unpolarized electron beams. However, with the experimental facility at the ILG, 
it will be possible to use polarized beams which will lead to an increase in the cross sections by a factor 
of four. 


Vs 

(GeV) 

(jpi (fb) 

(JP2 (fb) 

(JP3 (fb) 

crp4 (fb) 

Mb 

TC 

SG 

TC 

SG 

TC 

SG 

TC 

SC 

350 

0.067 

0.004 

0.048 

0.001 

~ 0 

~ 0 

0.132 

0.005 

5 

500 

0.453 

0.015 

0.399 

0.004 

0.006 

~ 0 

0.718 

0.025 

22 


TABLE III. Background cross sections for the signal process e“e“ —> W~W~ —>■ 4j after TG and SG for the 
four processes mentioned in the text for y/s = 350 GeV and 500 GeV. The number of background events A/s is 
computed at an integrated luminosity of 500 fb“^. 


For completeness, we also show the partonic cross sections for all possible decay modes of the IV-pair 
in Table IV. Due to larger branching ratios for the hadronic decay of W, the e“e“ —)■ 4j has the largest 
cross section. Note that for the numbers in Table IV, we have assumed the two VFs to be on-shell. The 
numbers for the second and fourth column will not change much if any of the IVs go off-shell; however, 
the numbers in the third column are valid only for on-shell Ws. 


Mn (GeV) 

4j (fb) 

2£- + $ (fb) 

r +2j + $ (fb) 

250 

0.075 

0.009 

0.055 

400 

0.084 

0.010 

0.059 


TABLE IV. Barton level signal cross sections of different channels arising from e“e“ — >■ W~W~ after trigger 
and selection cuts for illustrative cases of the heavy neutrino masses Mn = 250 and 400 GeV and with mixing 
parameter |I4jv| = 0.04 at y'i = 500 GeV. 


IV. MULTIVARIATE ANALYSIS 


In this section, we present a multivariate analysis for better signal-to-background discrimination, 
which leads to a better significance. For MVA, we use the Boosted Decision Tree (BDT) algorithm in 
the TMVA [110] framework. In the cut-based analysis for the channel e+e” —>■ Ni/^ —>■ e + 2j + $ 
presented in Section HI A, we have shown that two simple invariant mass cuts [cf. Eq. (15)] are sufficient 
to isolate the signal from the background with > 5cr significance, particularly for larger mixing (jVeTvj ^ 
0.02) and in smaller mass range (< 350 GeV). To probe smaller mixing (|I4Ar| ;< 0.02), we need to 
find an optimized set of cuts for the CBA in order to achieve 5cr significance. For this purpose, the 
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MVA technique is a useful tool to obtain the best sensitivity for a given set of parameters. For the 
MVA in the present case, we use nine variables, as listed in Table V with their relative importance in 
the MVA output response. These variables are chosen by comparing the signal trained for Mjv = 250 
GeV with |t4Af| = 0.04 at ^/s = 500 GeV with the background distributions. Each of these variables 
has some reasonable amount of discriminating power, as shown in Figure 18, where we have plotted the 
normalized signal and background distributions. It is important to mention here that this set of nine 
variables used here might not be the optimal one and there is always a scope to improve the analysis 
with cleverer choices of variables. In our analysis, we have used these simple kinematic variables which 
are less correlated and have sufficiently good discriminating power. We use the MVA input-variables set 
optimized for the parameters Mjq = 250 GeV, |I47v| = 0.04 and ^/s = 500 GeV for all parameter choices. 
In a more dedicated MVA, one can use different set of variables for different parameter points to further 
improve the analysis. 


Variable 

Importance 

Variable 

Importance 

Variable 

Importance 

PT{e) 

9.123 X 10"^ 

\riih)\ 

6.704 X 10“'' 

AR{e,ji) 

5.662 X lO""^ 

\vie)\ 

7.824 X 10"^ 


1.527 X 10"^ 

AA(ii,i2) 

6.180 X 10“'' 

PtUi) 

7.766 X 10"^ 

M{e,ji,j 2 ) 

2.522 X 10'^ 

Total $ 

1.625 X 10“^ 


TABLE V. Input variables used for MVA with Mn = 250 GeV, |14jv| = 0.04 and y/s = 500 GeV and their 
relative importance. This might vary for other sets of parameters. 



P.,(e) 



Input variable; T)(e) 
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linput variable: M(e,mE) 





Input variable: AR(e,j) 




linput variable: Total mE| 


LU 

i 0.07 


- 

Z 0.06 


- 

Z 0.05 



0.04 



0.03 


- 

0.02 


- 

0.01 
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FIG. 18. Signal (blue) and background (red) distributions of the input variables used for MVA. All the distribu¬ 
tions are drawn for Mn = 250 GeV, \VeN\ = 0.04 and y/s = 500 GeV. 

For MVA, we must always be very careful about overtraining the signal/background, since it can 
happen without the proper choices of the algorithm specific tuning parameters. Using the Kolmogorov- 
Smirnov (KS) test one can check whether a test sample is overtrained or not. Generally, the test sample 
is not overtrained if the KS probability lies within the range 0.1 to 0.9. For most cases, a critical KS 
probability value greater than 0.01 [111] ensures that the samples are not overtrained. In Figure 19, 
we have shown the KS probability values for the signal and the background of the BDT response and 
this confirms that both signal and background samples are not overtrained. We have ensured that we 
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TMVA overtraining check for ciassifier: BDT 





BDT response 


FIG. 19. BDT overtraining check for the parameter set Mm = 250 GeV, 14iv = 0.04 and = 500. 


do not encounter overtraining for the full parameter range. We have used two statistically independent 
samples for MVA, one for training and the other for testing. As can be seen from Figure 19, the signal 
and background samples are well separated in this BDT output, and hence, using an appropriated BDT 
cut, we can significantly increase the signal significance. 

In Table VI, we compute the significance for different Mm and fixed |14Af| = 0.04 at ^/s = 500 GeV 
with C = 100 fb“^. In Table VII, the significance is computed for different mixing values for a fixed 
Mm = 250 GeV at -y/s = 500 GeV with £ = 100 fb“^. From these two tables, it is evident that 
MVA is more powerful to probe smaller mixing and larger masses where the performance of the CBA is 
comparatively weaker. 


Mm 

^w/o cut 


BDT cut 

Ns 

Nb 

Ns 

(GeV) 

(fb) 

(fb) 

value 



VNs +Nb 

100 

39.28 

9.917 

0.162 

926 

301 

26.4 

150 

38.40 

22.61 

0.162 

1778 

946 

34.1 

200 

34.06 

22.73 

0.159 

1754 

786 

34.8 

250 

30.22 

21.94 

0.156 

1764 

651 

35.9 

300 

25.79 

19.87 

0.136 

1659 

475 

35.9 

350 

20.42 

16.33 

0.167 

1315 

431 

31.5 

400 

13.95 

11.43 

0.232 

876 

294 

25.6 

450 

6.605 

5.000 

0.218 

297 

194 

13.4 

475 

2.815 

2.361 

0.280 

119 

357 

5.5 

SM 

1022.5 

638.3 

- 

- 

- 

- 


TABLE VI. Partonic cross sections and cross sections after TC as defined in Section III A for 

the CBA with detector effects included for the process e''"e“ —>■ Nve e2j + $ a.t y/s = 500 GeV with mixing 
IVeivI = 0.04. Numbers of signal (Ns) and background (Nb) events for different Mm after passing the BDT cut 
are shown for luminosity £ = 100 fb“^. 

We also analyze the two subdominant channels e+e” —)■ Nv,. —>■ Zv^Vf. 2j + ^ and e+e” —>■ Ni/f, —> 
HVf^Vf. —26 + .^ for a few illustrative masses using a simplified MVA. For these two channels we have 
used eight variables, viz. priji), h(ji)|, tr, \'n{$T)\, M(ji,j 2 ), M($,ji), Ar]{$,ji) and scalar sum Hr 
of the pt’s of all visible particles, to distinguish signal from background. Before the training and testing 
of the samples, they were made to pass the trigger cuts given in Section III A. In Table VIII for the 
channel e+e“ —>■ Niy,, —?> ZveVf. —>■ 2j + we show estimates for significance for three illustrative mass 
values, viz. Mm = 100, 200 and 300 GeV for a fixed mixing |I4Af| = 0.04 at -y/s = 350 GeV. Similarly, in 
Table IX we show analogous results for the channel e+e” —)■ Nv^, —t Hve.Ve 2h + $ with mass values 
Mm = 150,200 and 300 GeV for a fixed mixing |Ve 7 v| = 0.04 at ^/s = 350 GeV. It is clear that using 
MVA, one can even probe these two subdominant channels with a good significance. It is important to 
note here that the significance for the process e+e“ —>■ Nv^, —>• 26 + .^ is subdued because of the 
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lUivI 

^w/o cut 

(fb) 

(fb) 

BDT cut 

value 

Ms 

Mb 

Ms 

V Ms -t Mb 

0.030 

17.00 

11.90 

0.158 

938 

580 

24.1 

0.025 

11.81 

8.27 

0.201 

582 

360 

19.0 

0.020 

7.55 

5.29 

0.227 

325 

00 

13.6 

0.015 

4.25 

2.98 

0.220 

187 

303 

8.4 

0.010 

1.89 

1.32 

0.336 

47 

60 

4.6 

0.0075 

1.06 

0.74 

0.309 

28 

71 

2.8 

0.0050 

0.47 

0.33 

0.359 

8 

30 

1.3 


TABLE VII. Same as Table VI for different mixing IVeivI and fixed Mjv = 250 GeV. 


smaller branching ratio of Hv than that oi N ^ Zv ior the benchmark points under consideration 
(see Figure 1). Also we consider a single 6-jet tagging in order to identify the Higgs, and this further 
reduces the signal due to the 6-tagging efficiency, which could vary from about 40% for low pT 6-jets to 
about 70% for high pt( 6). By demanding two 6-tags, the significance will decrease even further. Once 
the experiments start running and the detector performances are better understood, one can try to look 
for cleverer variables to improve this significance. 


Mm 

^w/o cut 


BDT cut 

Ms 

Mb 

Ms 

(GeV) 

(fb) 

(fb) 

value 



-h 

i 

100 

8.27 

5.89 

0.026 

O 

00 

8555 

5.1 

200 

13.61 

10.22 

0.017 

724 

6938 

8.3 

300 

4.39 

3.738 

0.082 

275 

5464 

3.6 

SM 

290.0 

218.9 

- 

- 

- 

- 


TABLE VIII. Same as Table VI for the process e''"e —>■ Nu^ —>■ 2j + $ at = 350 GeV. 


Mm 

^w/o cut 

^TC 

BDT cut 

Ms 

Mb 

Ms 

(GeV) 

(fb) 

(fb) 

value 



-h 

i 

150 

2.24 

0.795 

0.243 

36 

260 

2.1 

200 

5.49 

1.952 

-0.058 

149 

O 

00 

uo 

3.6 

300 

2.79 

1.045 

-0.005 

82 

1621 

2.0 

SM 

290.0 

44.45 

- 

- 


- 


TABLE IX. Same as Table VI for the process e^e —>■ —>■ HueVe 2h + $1 (only one 6 -jet is tagged among 

two 6 ’s) for c.m. energy ^/s = 350 GeV. The corresponding trigger cuts are defined in Section III A. 


V. CONCLUSION 

The seesaw mechanism provides a simple understanding of the observed smallness of the neutrino 
masses. Heavy Majorana neutrinos are an essential ingredient of the simplest seesaw scenario, and 
therefore, it is important to explore their experimental signatures in current and future experiments. In 
this work, we have studied the prospects of heavy neutrino searches at future lepton colliders. Due to 
their relatively cleaner environment, as compared to hadron colliders, lepton colliders provide us with a 
unique opportunity to probe the heavy neutrinos with higher sensitivity, which might shed light on the 
origin of neutrino mass and the associated new physics beyond the SM. In view of the increasingly hopeful 
case for an International Linear Collider in near future, we have mainly focused on various heavy neutrino 
production mechanisms at an e+e“ collider and have performed a detailed detect or-level simulation for 
the most promising channels. In particular, the production process e+e“ —>■ Nv leading to the final state 
e-|-2j-|-.^ has the largest cross section, and hence, we analyze this possibility in detail, taking into account 
realistic detector efficiency and SM background estimates. We have followed two different approaches 
to separate the signal from background. First, we present a conventional cut-based analysis using the 
distinct kinematic features of the signal over background. We show that in this approach, a light-heavy 
neutrino mixing parameter |I4Arp ~ 10“^ can be probed at 95% C.L. for a heavy neutrino mass up to 
400 GeV at ^/s = 500 GeV ILC with luminosity of 100 fb~^. Then we use a multivariate analysis with 
a boosted decision tree algorithm to make a better discrimination of the signal from background, which 
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enables us to probe even smaller mixing values up to IKatP ~ 10“^ or so for a larger range of heavy 
neutrino masses. 

The dominant production mode e“'"e“ —)■ Nu does not probe the Majorana nature of heavy neutrinos, 
which is a crucial aspect of the seesaw mechanism. To overcome this shortcoming at lepton colliders, we 
have studied a new production channel: e“'"e“ —)■ , leading to a lepton number violating same- 

sign dilepton plus four jet final state. Although the production cross section for this process is smaller 
than that oi e^e~ ^ Nv, the clean environment of e“'"e“ colliders makes this process almost background 
free. Thus, any positive detection in this channel would be strong evidence for lepton number violation, 
with profound implications for our understanding of the neutrino mass mechanism. We find that a 
detection of 10 signal events in this channel for \VeN\ = 0.04 would require an integrated luminosity of 
about 700 fb“^ for = 500 GeV, which is certainly within reach of next generation lepton colliders. 

Finally, we also study the lepton number violating process e“e“ —)■ W~W~, mediated by a t-channel 
heavy neutrino exchange, followed by the hadronic decay of the Ws. This channel is again almost 
background-free. Hence, a positive observation of this signal will be an unambiguous probe of the 
Majorana nature of the heavy neutrinos. We show that in this channel, we can achieve > 3a significance 
for the full range of the heavy neutrino mass considered here at a ^/s = 500 GeV ILG operating in the 
e“e“ mode with a luminosity oi C = 500 fb~^, provided the mixing is relatively large. With polarized 
beams, the signal cross section can be enhanced by a factor of four, thus leading to a better sensitivity. 

Although our numerical simulations have been done specifically with the ILG in mind, our analysis 
can be easily extended to other next generation lepton colliders, such as the CLIC, FCC-ee, GEPC or 
even muon-muon, electron-photon and photon-photon colliders, once the design parameters and detector 
specifics become clearer and readily available for these machines. 
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